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We present the results of the analysis of the AGN population in the deepest extragalactic hard X- 
ray survey. The survey is based on INTEGRAL observation of the 3C 273/Coma cluster region, 
and covers 2500 deg 2 with a 20-60 ke V flux limit 1 .5 times lower than other surveys at similar 
energies, resolving about 2.5% of the cosmic hard X-ray background. Using this survey, we 
can constrain in an unbiased way the distribution of hydrogen column absorption up to Nu= 
10 25 cm~ 2 . We put an upper limit of 24% to the fraction of Compton-thick objects. Compared 
to models of the AGN population selected in the 2-10 keV band, the LogrV-LogS diagram is 
generally in good agreement, but the Nh distribution is significantly different, with significantly 
less unabsorbed sources (Nn< 10 22 cm~ 2 ) at a given flux limit compared to the models. We also 
study the local hard X-ray luminosity function (LF), which is compatible with what is found in 
other recent hard X-ray surveys. The extrapolation of the 2-10 keV LF is lower than the hard 
X-ray LF. The discrepancy is resolved if AGN spectra typically present reflection humps with 
reflection fraction R ~ 1. Finally, we use the population properties of this survey to show that a 
future ultra-deep INTEGRAL extragalactic survey can result in a quite large AGN sample with 
enough objects at redshifts larger than z = 0.05 so that we can detect evolution in the hard X-ray 
LF. 
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1. Introduction 

The observation of the central parts of local galaxies have shown that a tight connection exists 
between supermassive black holes and their host galaxies [1, 2, 3]. Numerical simulations of 
galaxy mergers [4, 5] even suggest that the actively accreting supermassive black holes - what we 
call active galactic nuclei (AGN) - may play a fundamental role in the evolution of the galaxies. 
AGN surveys have therefore become an important part of observational cosmology. 

The X-ray domain is very efficient at targeting AGN, in large part because the vast majority of 
high-latitude X-ray sources are AGN. However, it is expected that X-ray surveys may provide only 
a partial view of the AGN population. The HEAO- 1 satellite discovered an apparently diffuse X-ray 
emission at high galactic latitude [6] whose spectrum peaks around 30 keV. While AGN could be 
the most important contributor to the so-called cosmic X-ray background, the spectrum could not 
be explained by the sum of AGN spectra, unless there exists a population of highly absorbed AGN, 
with hydrogen column densities Nu larger than 10 22 cm~ 2 , and even Compton-thick objects with 
jVh> 10 cm~ 2 [7, 8, 9, 10]. These Compton-thick objects emit very little radiation below lOkeV 
and thus require deep X-ray observations. In a recent detailed modeling of the AGN population 
based on the known AGN population up to z ~ 3 and its extrapolation to higher redshifts, [1 1] found 
that the population of Compton-thick AGN should be as large as that of moderately absorbed AGN. 

Absorption is much less efficient in the hard X-ray domain (> 20keV) than in lower X-ray 
bands. Surveys in the hard X-ray have therefore the potential of detecting bright AGN with minimal 
bias in the A^h distribution below ~ 10 25 cm~ 2 . INTEGRAL and SWIFT are two satellites with such 
survey capabilities which allowed to study the AGN population over the full sky [12, 13]. However, 
the sensitivity limit of these surveys remains very low compared to those of surveys conducted by 
XMM-Newton or Chandra. 

We present here the analysis of new INTEGRAL data covering a ~ 2500 deg 2 region of the sky 
centered around 3C 273 and the Coma cluster. Thanks to these data, this survey is the deepest high- 
latitude survey in the hard X-rays. Most results have been presented in detail in [14]. We report 
on the population properties of the detected AGN, and in particular their absorption properties 
compared to what is known from the AGN population observed in 2-10keV X-ray surveys. We 
show here what can be derived from the comparison of the LFs in the 2-10 keV and 20-60 keV 
ranges. Finally, we explore the usefulness of performing a much deeper INTEGRAL high-latitude 
survey in the form of a 5x5 dithering pattern with up to 20 Ms total exposure time. 

2. The deepest INTEGRAL/IBIS extragalactic field 

Several INTEGRAL core-programme and open-time observations have covered the sky region 
around Coma cluster and 3C 273, which is illustrated on Fig. 1. We selected all available INTE- 
GRAL pointings within 30 degrees of a position located between these two sources, which resulted 
in 1660 pointings for a total elapsed observing time of 3 936 234 s and a dead-time corrected good 
exposure of 2 733 202 s. Most pointings belong to four 5x5 dithering patterns repeated several 
times, plus a specific rectangular pattern used during the core-programme observation. 

Sky images in the 20-60 keV energy ranges were created from the data taken by the ISGRI 
detector of the IBIS imager on board INTEGRAL [15]. The 3000x3000-pixel resulting mosaic im- 
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Figure 1: INTEGRAL/IBIS all-sky survey [12] showing the region of the deepest high-latitude INTEGRAL 
exposure analyzed in [14] (blue area with red contours). 
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Figure 2: Properties of the mosaic of the 3C 273/Coma region. Left: Histogram of the pixels' significance. 
The red curve is a Gaussian fit to the part of the histogram with a < 3. Right: Surface of the sky over which 
a given flux results in a significance (7 > 5.5 (solid line) and a > 5.0 (dashed line). The vertical dotted line 
shows the sensitivity limit of other hard X-ray surveys. 
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age was built in equatorial coordinates with a tangential projection using a factor-2 over-sampling 
when compared to the individual input sky images; this results in a pixel size of 2.4 arcmin in the 
center of the mosaic and of about 1.6 arcmin in the outskirts of the image, roughly 40 degrees away 
from the center. 

To avoid any systematic bias in the pixel significance distribution, we excluded the area for 
which the effective exposure time turned out to be smaller than 10 ks. The resulting pixel signifi- 
cance distribution is shown on Fig. 2 left. The part of the histogram with significance smaller than 
3 is perfectly modeled with a Gaussian distribution. The sigma of the Gaussian is however 1.1, 
instead of the expected value of 1. The centroid is found at a value of 0.009, which has negligible 
impact on the analysis. 

We extract the fluxes and significance of the candidate sources by fitting a PSF directly over 
the mosaic. In order to derive the sensitivity as a function of effective exposure time, we extracted 
fake sources at random over the mosaic and estimated their significance distribution. Taking into 
account repetitions (our mosaic contains about 13 800 independent pseudo-sources), we obtain that 
a source with a 5a significance has only 20% chances of being real, and the figure increases to 
85% if we set the threshold to 5.5a. 

3. Properties of the hard X-ray AGN population 

We find 34 candidate sources in the mosaic with significance larger than 5a. Source param- 
eters are listed in [14]. Taking into account the probabilities that some sources are fake, we end 
up with approximately 26 real sources. 22 of these sources are found to coincide with known 
low-redshift (z < 0.05) AGN. We found that the probability of such chance coincidence is of the 
order of 1%. Thus we can consider that these 22 sources are real. We nevertheless use the formal 
probabilities of the sources being real in all analyzes of the AGN population below. 

3.1 Hydrogen column density distribution 

17 of the 34 candidate sources have adequate X-ray measurements allowing us to determine 
their hydrogen column density, Ah- For the 17 remaining objects, we check whether there is a 
counterpart in the ROSAT all-sky survey bright source catalog (RASS-BSC) [16]. This catalog 
is quasi-complete over the full sky down to a flux of 0.05 ctss -1 in the 0.1-2.4 keV band. We 
calculate the expected count rate in the ROSAT band for different values of intrinsic Ah assuming 
a power-law intrinsic emission with index T = 1 .9. While the method is rather crude, it is sufficient 
to obtain a moderately accurate estimate of Ah- In particular, we checked that it is able to identify 
correctly, for all 17 sources with A^h measurements, whether an AGN is absorbed (Nu> 10 22 cm -2 ) 
or not. When the source is not detected in RASS-BSC, A^h is roughly constrained within the range 
2 x 10 22 -10 25 cm" 2 . 

Figure 3 shows the distribution of the intrinsic hydrogen column density Nu for the 34 sources 
detected in the mosaic and weighted by their probabilities. The fraction of absorbed objects, i.e. 
those with Nu> 10 22 cm~ 2 is found to be 70%; if one discards the 1 1 sources without X-ray coun- 
terparts, this figure becomes 46%, making it a stringent lower limit. None of the 23 sources with 
measured or estimated Nu are Compton-thick (A f H> 10 24 cm -2 ), placing an upper limit to the frac- 
tion of Compton-thick objects of 24%. These figures are consistent with those found in previous 
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Figure 3: Histogram of the A^h distribution for the candidate sources. The blue histogram shows sources 
with measured intrinsic N\\S. The green histogram shows the sources whose is derived from the presence 
of a counterpart in the RASS-BSC. The red histogram shows the sources without counterparts in the RASS- 
BSC, with log A^h evenly distributed between 22.25 and 25. 



hard X-ray surveys. In a recent modeling of the cosmic X-ray background [11], it was predicted 
that, at the level of 10~ n erg s _1 cm~ 2 , the fraction of absorbed AGN (with Nn> 10 22 cm~ 2 ) 
should be 65%, which is perfectly compatible with our measurement. Their expected fraction of 
Compton-thick AGN (15%) is also compatible with our upper limit, although the lack of detection 
of any true Compton-thick object makes our constraint rather weak. 

We do not observe any drop in the fraction of absorbed objects at high luminosity, contrarily 
to several observations in the 2-10 keV range [17, 18, 19]. While this may due to the small hard 
X-ray sample sizes in general and to the very small number of high-luminosity objects in particular, 
we point out that currently no hard X-ray survey has observed this effect with some significance. 

3.2 Comparison with the AGN population observed at 2-10 keV 

In Fig. 4 top left we compare our source counts to the predictions of the 2-10 keV population 
model of [17]. We predict the 20-60 keV source counts by integrating the 2-10keV model AGN 
population over the < z < 2, 10 41 < L2-i0keV < l^ 47 er § s range. The conversion from rest- 
frame intrinsic (i.e. before absorption) 2-10 keV luminosity to observed frame 20-60 keV flux is 
made using the same spectral model as in [17], namely a power-law spectrum with T = 1.9, a cut- 
off rest-frame energy of 500 keV, and a reflection component from cold material. With this spectral 
model, S2o-60keV = 1-24 x L2-iokev/47T<iL at redshift ~ (Jl is the luminosity distance in cm). 
We use here the model of [17] which includes a mix of unabsorbed and Compton-thin sources, 
but no Compton-thick sources; the Compton-thick population has indeed not been measured and 
has been treated somewhat arbitrarily by adding a number of these objects equivalent to that of the 
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Figure 4: logA^-logS diagrams for the sources detected by INTEGRAL. In all panels the solid lines are 
from our work and the dashed or dotted line are model extrapolations from [17]. Top left: logA^-logS for 
the full list of candidate sources with the la uncertainties. The solid line is the power-law best fit. Top right: 
logAMogS in two redshift bins; red: z < 0.05; blue: z > 0.05. Bottom left: logA^-logS in two luminosity 
bins; red: L20-60keV < 10 43 ' 5 ; blue: L20-60keV > 10 43 5 . Bottom right: logA^-logS in two Na bins; red: 
N H < 10 22 cm- 2 ; blue: N H > 10 22 cm- 2 .TJ03' refers to [17]. 



Compton-thin ones. We have therefore assumed that absorption effects on this population in the 
20-60 keV band are negligible. 

Under these assumptions, we can see that the 2-10 keV model provides a good match to both 
the slope and normalization of the total 20-60 keV source counts (see Fig. 4 top left), leaving little 
room for a significant additional population of moderately Compton-thick sources. We examine 
the population in more detail by splitting the sample into low (z < 0.05) and high redshift sources, 
low (L2o-60keV < 10 43 5 erg s _1 ) and high luminosity sources, and low (Nu< 10 22 cm~ 2 ) and high 
absorption sources. In Fig. 4 top right we show the source counts and model predictions for red- 
shifts either below, or above 0.05. Redshifts are unknown for 6 sources. The observed counts and 
model predictions agree reasonably well given the relatively small numbers of observed sources. 
Above 52o-60keV ~ 10~ n erg cm~ 2 s _1 , the 2-10keV model from [17] overpredicts the number of 
observed sources with z > 0.05. 

In Fig. 4 bottom left we show the source counts and model predictions for observed sources 
above and below a luminosity of L2u-60keV = 10 43 ' 5 erg s _1 . For our given spectral model this 
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Figure 5: Hard AGN X-ray luminosity function from the 3C 273/Coma survey. The black points use the 
V /Vmax estimator, while the black curve has been obtained with a maximum likelihood. The blue curve is 
the LF from an INTEGRAL all-sky survey [12] and the red curve is from the Swift/BAT survey [13]. The 
LF from the luminosity-dependent density evolution model found in the X-rays [17] and converted to the 
hard X-ray domain using T — 1 .9 is shown in green. 

corresponds to L2-iokeV = 10 43 ' 41 erg s _1 (as described above). The observed counts and model 
predictions are roughly in agreement, given the small number of sources. 

Absorbing column estimates or lower limits are available for our entire sample and so, in 
Fig. 4 bottom right, we show the source counts and model predictions separated into sources with 
Nu either greater than, or less than 10 22 cm 2 . Here it is clear that the 2-10 keV model is a poor 
predictor of the observed source counts. We see that the more absorbed sources constitute at least 
2/3 of the total 20-60 keV source counts over the flux range of the sample, whereas the 2-10 keV 
model of [17] predicts equal number of absorbed and non-absorbed sources over the luminosity 
and redshift range probed by the INTEGRAL observations. Taken at face value, it means that the 
observed and model Nr distributions differ significantly. A follow-up with more sensitive medium 
energy X-ray observations is therefore needed to confirm this result, especially for the sources with 
no direct A^h measurement. 

4. The hard X-ray luminosity function 

We use the candidate sources firmly detected as AGN to determine the AGN luminosity func- 
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tion (LF) in the 20-60 keV energy range. We use only objects with redshifts < z < 0.05 in order to 
study the local AGN population with negligible evolution. Because of the small number of sources, 
we use a parametric method to derive the LF. We assume the standard AGN luminosity function: 

*m= ( 4i) 

dlogL (L/L*) a + (L/L*)P' 

which describes a broken power-law, changing from index a to index /3 at characteristic luminosity 
L* . The parameters L* , a and /3 parameters are determined using a maximum-likelihood (ML) test, 
based on the idea proposed by [20], modified to take into account our sensitivity map. ML-based 
methods are insensitive to sample sizes, therefore <J>* cannot be estimated this way. Thus we also 
estimate the LF using the standard non-parametric V /V max estimator. <£* is then obtained by fitting 
the V /V max LF estimates with the <&(L) distribution from Eq. (4.1), letting only <I>* free and fixing 
the other parameters to the ML values. We obtain the following results for the four LF parameters: 

logL* = 43.66^ ergs- 1 
oc = 0.85 ?« 

(4 2) 

J8 = 3.12^ (4 - 2) 
<£* = 1.121^ 10- 5 Mpc- 3 

As can be seen from Fig. 5, the AGN hard X-ray luminosity functions from INTEGRAL all-sky 
survey [12] and the Swift/BAT survey [13] are perfectly compatible with our estimate. 

<£* and L* are very strongly correlated; therefore the LF normalization is known with a much 
better accuracy than the uncertainty on <I>* suggests. As an example, the luminosity density inte- 
grated above L 2 o-60keV = 10 41 erg s _1 is: 

W(L>10 41 ) = /-L 5 g I dlogL= 

= 0.90+°;^ 10 39 erg s" 1 Mpc~ 3 

The luminosity density is again perfectly compatible with [12] and [13] (W(L > 10 41 ) = 1.02 10 39 
and 1.03 10 39 erg s~ ! Mpc~ 3 respectively). 

Figure 5 left also shows the local, z = 0, X-ray LF from [17]. This LF is significantly below 
the hard X-ray ones by a factor 1.5. There are two extreme ways to reconcile the LFs: either the 
X-ray LF is shifted vertically, which correspond to a change in AGN density, or horizontally, which 
corresponds to a change in AGN luminosity. 

A change in density can be explained if there is a fraction of objects that have been completely 
missed in X-ray surveys, for instance the Compton-thick population. This would imply that one 
third of the AGN in our sample are Compton-thick, which is above our upper limit. 

Alternatively, a change in luminosity resulting in a horizontal shift in the LF might be ex- 
plained if we underestimated the hard X-ray luminosity of AGN observed in the X-rays. This is 
perfectly possible, since we had to transform the unabsorbed 2-10 keV X-ray luminosity assum- 
ing a spectral shape, which we chose as a simple power-law with photon index T = 1.9. The LFs 
would become consistent if we choose instead a photon index T = 1.7. Interestingly, we obtain a 
very good approximation to a power-law with photon index T = 1.7 if we combine a power-law 
with photon index T = 1 .9 with a reflection with reflection fraction R = 1 . Thus the comparison of 
the LFs is evidence that large reflection fractions with R ~ 1 are common in AGN. 
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Figure 6: Surface of the sky over which a given flux results in a significant detection. The red curve 
corresponds to a 5x5 dithering pattern with a 10 Ms total exposure time. The blue curve corresponds to the 
same observing pattern, but with a 20 Ms exposure time. The black curve is the 3C273/Coma survey [14] 
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Figure 7: Simulated LogiV-LogS diagram (red) resulting from a 5x5 dithering pattern with a 10 Ms total 
exposure time (left) and a 20 Ms total exposure time (left) respectively. The blue line show the 1 a uncer- 
tainties. The green line is a power-law fit to the Log YV-Log S diagram. 



5. Prospects for an ultra-deep INTEGRAL extragalactic survey 

In terms of sensitivity, coded-mask telescope performances are well below those of focusing 
telescopes. While telescopes able to focus hard X-rays will fly in a few years and will make 
wonderful contributions to the study of the AGN population, their fields of view are much smaller 
than those of Swift/BAT or INTEGRAL/IBIS, which limits their survey capabilities. It is therefore 
worth exploring the prospects of an ultra-deep INTEGRAL/IBIS extragalactic observation that 
could realistically be completed before the end of the mission. 
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We simulated INTEGRAL/IBIS observations following a 5x5 dithering pattern with 10 Ms and 
20 Ms total exposure times. Fig. 6 shows the surface of the sky over which a given flux sensitivity 
would be reached with such observations. On the basis of these curves, we can draw a population of 
AGN assuming that the Log A^-Log 5 diagram obtained in the 3C 273/Coma survey. Log Af-Log 5 
diagrams of such simulated populations are shown on Fig. 7. 

The average number of sources is 57 for the 10 Ms observation and 95 for the 20 Ms one, with 
0-area minimum flux limits of 3 and 2 x 10 erg s _1 cm~ 2 respectively between 20 and 60keV. 
Below 10~ n erg s _1 cm~ 2 , AGN with redshifts z > 0.05 dominate the full population. We find 
that at least 20 objects with z > 0.05 will be found in the 10 Ms observation, and at least 50 in the 
20 Ms one. 

We can translate these sample sizes in uncertainties on the normalization of the luminosity 
function. We obtain an uncertainty better than 22% with the 10 Ms observation and better than 
14% with the 20 Ms one. The X-ray luminosity function normalization increases very rapidly, like 
(1 +z) 4 , at very small redshifts [17], implying a 20% density evolution between z = 0.025 and 
z = 0.075. As the local hard X-ray luminosity function will be quite well constrained at the end 
of the Swift mission, this means that an ultra deep INTEGRAL/IBIS extragalactic observation has 
the potential to constrain AGN density evolution in the hard X-rays. 

6. Conclusion 

While the extragalactic hard X-ray survey presented here covers an area much smaller than 
other recent similar surveys [12, 13], we reached limiting fluxes a factor 1.5 times deeper. We 
obtained a source density of 0.013 deg" 1 at a 20-60 keV flux level of 10 -11 erg cm s . We 
resolved approximately 2.5% of the cosmic X-ray background. Comparison with the 2-10 keV 
domain shows that the total LogA^-LogS diagrams are compatible, indicating that there isn't any 
large population of bright Compton-thick objects missed in the 2-10 keV surveys and appearing in 
the hard X-rays. 

No object is our sample is truly Compton-thick, although the upper limit to the Compton-thick 
fraction is compatible with the fraction expected from models of the cosmic X-ray background [11]. 
Nevertheless the absence of such objects makes that the case for the existence of such population 
at the required level is rather weak. Follow-up observations of our sources without adequate Ah 
measurements are in progress, and may solve the puzzle. It seems however quite probable that we 
shall end-up with conflicting Ah distributions, which will have to be explained. 

We present a truly local hard X-ray luminosity function (z < 0.05). We find a LF quite com- 
patible with the latest all-sky INTEGRAL and SWIFT surveys. A discrepancy is however found 
between the hard X-ray LFs and the 2-10 keV ones, which we interpret as an excess luminosity in 
the hard X-rays, which could result from a reflection hump with reflection fraction R ~ 1 . 

While still suffering from low sensitivities compared to their modern counterparts working in 
the soft and medium X-rays, INTEGRAL and Swift are the only instruments currently available 
to perform surveys above 15keV. The importance of this energy domain, which is unaffected by 
obscuration below ~ 10 25 cm~ 2 , is such that efforts to build statistically representative samples 
must be pursued. While the Swift mission will cover the full sky with quite deep observations, it is 
possible for INTEGRAL to perform ultra-deep observations way beyond what could be obtained 
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with Swift. We showed here that such INTEGRAL observations would produce numerous faint 
sources allowing us to quantify the AGN population with very good accuracy. In addition they 
would allow us to explore the evolution of the AGN population in the hard X-rays for the first time. 
Therefore ultra-deep INTEGRAL high-latitude observations, maybe reaching some 20 Ms total 
exposure time, have a very strong scientific potential and should absolutely be performed before 
the end of the INTEGRAL mission. 
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